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^ . Abstract 

Determining the spin of new particles is critical in identifying the true theory among various 
extensions of the Standard Model at the next generation of colliders. Quantum interference 
between different helicity amplitudes was shown to be effective when the final state is fully 
reconstructible. However, many interesting new physics processes allow only for partial recon- 
struction. In this paper, we show how the interference effect can be unambiguously extracted 
even in processes that have two-fold ambiguity, by considering the correlation between two de- 
cay planes in e"'"e~ collisions. 

PACS numbers: 12.60.-i, 12.60.Jv, 14.80.Ly 
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The Large Hadron Collider (LHC) will soon usher us into the arena of the electroweak symme- 
try breaking scale and beyond. Once the TeV regime is explored, it is highly anticipated that the 
true theory for the origin and stability of the electroweak scale [1] will be revealed. 

One possible result from the LHC, generically predicted in new most models, is the presence 
of new particles partnered with some or all of the SM particles. For instance, every SM particle 
in the minimal supersymmetric standard model (MSSM) [21 [3] has a heavier partner whose spin 
differs by 1/2. Alternatively, in the minimal universal extra dimension (UED) model [4] with a 
compactified electroweak-scale extra space dimension, each SM particle is paired with a tower of 
Kaluza-Klein (KK) states with identical spin. Thus, model-independent spin measurements are 
crucial in discriminating among many extensions of the SM. 

There have been several proposals for measuring spin at both the LHC and the prospective 
e^e~ International Linear Collider (ILC) [5]. Threshold scans in e~^e~ collisions can be used to 
distinguish scalars from spinors at the ILC, as the scalar production cross section increases slowly 
~ f3'^ while the spinor cross section increases steeply ~ /3 [6]. [Such a method cannot be used at the 
LHC as the center of mass (cm.) energy at the parton level is not fixed.] The production angle 
can give insight on spin as well. The polar-angle distribution of a s-channel pair produced scalars 
is proportional to sin^ B, while for spinors it approaches 1 + cos^ O asymptotically at high energies. 
The presence of i/u-channel exchanges may render the production-angle measurement of spin more 
demanding [6] , although it is feasible in some cases [71 [8] . The polar-angle dependence in decays 
at the ILC [8] and the invariant-mass distributions in sufficiently long decay chains at the LHC [9] 
can also be used for spin measurements. However, these techniques rely strongly on the final state 
spins and the chiral structure of couplings. 

In this report we study the fully- correlated azimuthal-angle distributions in the production of a 
new particle-antiparticle pair in e~^e~ collisions and both of their sequential decays |10j . [This work 
is a natural extension to the previous works [11^ [T2] where the azimuthal-angle correlation of the 
production and only one of the decays has been investigated.] These distributions develop through 
quantum interference between the different helicity states in a coherent sum. By extracting this 
angular dependence, we can determine which helicity states contribute to the sum, and thus the 
spin of the decaying particle in a model-independent way. To be specific, we restrict ourselves to 
production of an electrically-charged particle-antiparticle {F'^F~) pair in e~^e~ collisions and decay 
of each produced particle to a charged particle and an invisible particle X; 

e+e- ^ F+F- ^ (/+x)(/"x) - ^. (1) 

As suggested by the WIMP solution to the cold dark matter puzzle, large missing energy signatures 
are considered likely at the TeV scale and are generic in many extensions to the SM|j 

For example, the lightest supersymmetric particle (LSP), typically the lightest neutralino Xii in 
supersymmetric (SUSY) models with i2-parity; and the lightest KK odd particle (LKP), typically 
the lightest KK gauge boson 71, in UED models with KK parity are stable and escape detection. 
Prototype processes in the SUSY and UED models with the same event topologies as the process ([1]) 

*The SM processes e"'"e~ — > t^t~ — > (7r+i>T-)(7r~^'T-) and e^e~ — > W^W~ (t^ vt)(t.~ vi) with t = e and /i also 
have a large missing energy signature carried away by invisible neutrinos. 
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are 

e+e- ^ /i+ /i- ^ (/^+X?)(/^-X?) ^ ^, (2) 

e+e" ^ /i^i/ifli ^ (/i+7i)(^"7i) ^ V- (3) 

Both generate the same experimental signatures [i^ ijT Jfi with large missing energy. 




Figure 1: The definition of the polar angles of the visible particle f momentum in the rest 
frame of the decaying particle and of the correlated azimuthal angle (j) between two decay planes 
formed in the correlated production- decay process X F~F~^ — > {f~x)if~^x) ^'^ '"est frame of 
the X = {e'^e~} system, corresponding to the e~^e~ cm. frame for the processes considered in this 
report. Here, X = {e~^e~} denotes any single- or multiple-particle intermediate state formed in 
e+e~ annihilation. Note that (p is invariant under the Lorentz boost along the F^ flight direction. 

The characteristic observables for measuring spin of the particles F^ through the process ([T]) 
are the angular distributions of the final-state particles in the F^ decays, encoding the helic- 
ities of the F^ states. We denote the polar angles of the particles in the rest frames of the 
F^ particles by 9^, and the azimuthal angles by with respect to the production plane defined 
by the e~ and F~ momentum directions, respectively. Then the angle (j) with its range [0, 2tt] 
between the two decay planes (see Fig.[T]) is the azimuthal angle defined by the angle difference 
(/) = <p*^ — (mod 27r) invariant under any Lorentz boost along the F^ flight direction. 

If we label the F^ helicities by X± and Aj_, the joint production-decay distribution reads: 

w{E,^;Q;eiA*±) = E v',r'/ {E,^,e)v~ (r (e;,<A;), (4) 

where Ecm is the e~^e~ cm. energy and is the production angle of F~ with respect to the e~ 
direction, and j is the spin of the particle F^. The production density matrix V is defined in terms 
of the helicity amplitudes T of the process e~^e~ — > F^F~ for unpolarized beams by 

cr±=±l/2 

where a± is the helicity, and each F^ decay density matrix T)^ has a simple azimuthal-angle 
dependence of a pure kinematical origin as 

^A^±A^(^i,0i) = ^tA'J^i) e^'(^±-^±)^i, (6) 
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reflecting an overall rigid rotation of the decay plane around the parent particle momentum. 



Integrating the joint production-decay distribution W in Eq. Q over the production angle 
@, the decay angles and (j)'^ with the azimuthal angle (j) fixed, we can derive the correlated 
azimuthal-angle distribution between the two decay planes as 



dC 



d cos d cos 6*_ d cos 9* 



(7) 



We note from Eqs. ([4]) and ([6]) that the dependence on (j)*^ is of the form exp[— i(/)!j_(A_— A'_— A++A'_,_)] 
so that the integral over (j)*^ leaves only those terms in Eq. Q satisfying the relation A+ — A^ = 
A_ — A'^ = A in the range [— 2j, 2j]. If the distribution is further integrated over the angle cj)^ only 
the incoherent terms with \± = Aj_ survive. Thus, any non-trivial azimuthal-angle distribution 
indicates the presence of quantum interference between the different helicity amplitudes. 



In weakly-interacting and CP-invariant theories with negligible particle-width and loop effects, 
the general form of the normalized azimuthal-angle correlation for the production and decay of a 
spin-j particle paiiQ is 



IdC 
CM 



27r 



[1 + Ai cos(0) + 



A2j cos{2j(j))] . 



(8) 



Each coefficient can be worked out from the standard rules of constructing matrix elements. How- 
ever, it is guaranteed on a general footing [S] that the highest non-vanishing coefficient is always 
A2j. This is because the production of a charged pair in e+e~ collisions gets a non-zero spin-1 
photon-exchange contribution to the production amplitudes, %T_a+;±j±j- However this term tends 
to be suppressed by ~ m?p±/E'^^ at high energies because of a final-state helicity flip. Thus, the 
spin j can be determined by identifying the highest cos{2j(p) mode at a cm. energy not far away 
from the production threshold, if the production amplitudes contributing to the coefficient A2j are 
not so suppressed and the decays —>■ f^x do not have too small polarization analyzing powers. 



The correlated azimuthal angle distribution for the smuon-pair process ([2]) is flat, as it must 
be, and the distribution for the KK muon-pair process ([3|) is given by 



1 dC 

Cd6 



1 

2^ 



1 



2m2^ 



cm + , 



1 -I- 2m2 /m^ 1 



COS ( 



(9) 



It is apparent from the expression ([9]) that (a) the coefficient of the highest cos 4> mode is maximal 
in magnitude at the production threshold and it decreases rapidly with increasing energy in con- 
formity to the general rule as outlined above, and (b) it is very sensitive to the values of the fiRi 
and 7i masses, leading to the restriction that the magnitude of the coefficient Ai of the highest 
cos(j) mode cannot be larger than 7r^/48 ~ 0.206. 



As well known [71 [8] , there exists a two- fold discrete ambiguity in completely reconstructing the 
F^ four-momenta and thus the azimuthal angles 4'± in the process ([1]) in the laboratory frame. 



^Even in the CP-noninvariant case, all the sine terms are washed out by taking the average over two possible 
azimuthal angles, which is unavoidable due to a two-fold ambiguity in reconstructing the F"^ momentum as described 
in the following. 
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Figure 2: An illustration for the relation (pp = 271 — (j)T between the true and false azimuthal angle 
differences, (pT = ^.''^d cpp, leading to the equality: coscpp = coscpx = cos</>. The dashed-line 
circles are the circles projected on a sphere of a unit radius (centered on the e^e~ interaction) 
of two cones satisfying the relation ilU\) . The solid black dots indicate two solutions for the 
direction and the unit vector n± stands for the direction. 

even if all particle masses are known. In [TT] it was shown that this ambiguity could obscure the 
helicity information in (p^, curtailing their use as measurements of spin. Nevertheless, the cosine 
of the azimuthal angle cp is unambiguously determined by measuring the four-momenta event 
by event. To prove this important point analytically, let the pair produced particles and the 
invisible particle x have mass m± and ttt-q and denote the flight direction in the laboratory 
frame by a unit vector n±, respectively. Then, the opening angles 6± between the visible tracks 
and the parent F^ particles in the laboratory frame can be determined from the relation 

ml-ml = E,^Ef± (^1- cos , (10) 

defining two cones about the f~^ and /~ axes which intersect in two lines - the true F^ flight 
direction and a false direction. True and false solutions are mirrored on the plane spanned by the 
and /~ flight directions, leading to the relation (px = 2Tr — (pp = cp between the true and false 
values, (pT and (pp, of the azimuthal angle (p (see Fig. [2]). Therefore, the cosine of the azimuthal 
angle is uniquely determined and its expression is given by the simple expression: 

cos (p = (n+ • n_ + cos 6^ cos 9^ ) / sin 6^ sin ^_ , (11) 

expressed in terms of the unit vectors, n±, and the opening angles, 6±. 

Two cosines, cos{na(p) and cos(nb4>), foi' &iiy integers and nj,, are functions of cos (p and orthog- 
onal to each other for Ua ^ rij,. Therefore, we can project out all the coefficients Ai^ (k = 1, . . . , 2j) 
by fitting the expression of (l/C) dC/dcp in Eq. ^ to the distribution measured experimentally. 

For a numerical demonstration of this spin-determination method, we compare the correlated 
azimuthal-angle (p distribution of the SUSY process ([2]) with that of the UED process ([3]) in a 
specific scenario, which we will simply call "BCMM" for convenience in the following, with the 
particle mass spectrum, 

BCMM : m± = m-± = m ± = 200 GeV and mo = m^o = m-y, = 50 GeV. (12) 
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We stress that the mass spectrum is chosen only as a simple illustrative example for SUSY and 
UED models with different spins but similar final states and so the spin-determination method 
demonstrated here can, in principle, be exploited equally for any other scenarios beyond as well as 
within the SM. 




Figure 3: (a) The cm. energy dependence of the production cross sections for smuons and KK 
muons; and (h) the normalized azimuthal- angle distribution {2t: /C)dC / d(j) in the BCMM scenario 
with Ecm = 450 GeV. The solid (dashed) line is the distribution without (with) the rapidity cut on 
the fj,^ directions and the total missing momentum, \ri\ < 2.5. 

The total cm. energy at the ILC is expected to reach up to 1 TeV, and an integrated luminosity 
of 500 fb~^ is not unrealistic. For the mass spectra chosen, we expect several thousand to several 
hundreds of thousands of events available as shown in Fig.[3||a). [At high energies we note that the 
fi^i scalar cross section scales in the same way as the jl^ spinor cross section, but with a coefficient 
4 times as large due to differences in the number of spin degrees of freedom, as familiar from QED 
processes.] 

To simulate the effects of experimental cuts which are unavoidable due to the geometry of the 
detector, we place cuts on the pseudo-rapidity of the fi^ directions and the total missing momen- 
tum: |?/| < 2.5, as otherwise the leptons would vanish unseen down the beam and the missing 
momentum is not guaranteed to be carried away by the invisible x particles. Two representative 
distributions for scalar and spinors (both with and without rapidity cuts) are shown in Fig. ^h). 
The SUSY distribution is fiat and the UED distribution has a clear cos (j) dependence above a flat 
distribution, as expected. Furthermore, the rapidity cut reduces the total number of events by 
about 2% but hardly modifies the correlated azimuthal-angle distribution. It is therefore apparent 
even at this level of analysis that the non-trivial correlated azimuthal-angle distribution contains 
the spin-1/2 information of the KK muon, /U^^. 

Using the least-square method we fit the generated distributions to (l + Ai cos(j) + A2 cos2(/))/27r 
after placing the cut on the pseudo-rapidities of the £^ and the total missing momentum. Only 
the coefficient Ai for the scalar flu and spinor fiRi are shown in Fig. HI because the coefficient 
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A2 is found to be extremely small (< 0.1%). This is consistent with the statistical errors intro- 
duced from the event generation and fitting procedure and confirms our general argument that 
Ai = for i > 1. As can be seen, in the BCMM spectrum the values of Ai for the scalar muons, 
fiR, are consistent with zero for all energies, as expected. For the KK muons, fi^^, the coefficient 
Ai is manifestly non-zero, allowing us to clearly distinguish the spinor field from spin-0 scalar states. 
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Figure 4: Coefficient Ai for the BCMM point as a function of the e'^e" cm. energy for both 
the SUSY fb^ and UED fi^^ pair production with an integrated luminosity of 500 fh~^ . Error bars 

f^^Xi) — -^^(/^Ri ~^ A'^Ti) = 1 correspond to the 1-a uncertainty range. 



obtained with Br{fi^ 



Due to the large Ai signal for the ii^^ process on the order of 10% and the negligible rapidity-cut 
effect on the normalized azimuthal-angle distribution, the ILC with a large integrated luminosity 
can discern that the particle is not a spin-0 but a spin-1/2 particle. However it is not always 
guaranteed that, in the case of higher spins, the highest mode is sufficiently large to be discerned. 
As a result, the question still remains whether this spin-determination method may be practically 
applied to discriminate spinors from vectors and vice versa in general cases. It is worthwhile to 
study the case of pair production of massive spin-1/2 fermions in SUSY contrasted with vector 
production in UED. 



In particular, we consider the processes - the production and two-body decays of a lighter 
chargino pair and those of a first KK VF-boson pair: 

(£+p,)(rp;), (13) 



e ' e 



Xi Xi 



e+e- ^ W+W{ ^ (^+z^£i)(r%), (14) 

where the charged leptons can be of either a muon or an electron type. An explicit calcula- 
tion shows that, when CP is preserved and all the absorptive parts are negligible, the correlated 
azimuthal-angle distributions (27r/C) dC/d(j) are given in terms of the production matrices by 

^2 



Xi Xi 



1 + 



TT 



^e{pZZ) cos(0). 



(15) 
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32 \2m^^±+m^^ 

2 2 \ 2 



3?e(plt + cos(20), (16) 



where the super /sub-scripts =b stand for the hehcities ±1/2 and ±1 of xf and W^, respectively. 
The integrated production density matrix p is defined as|i| 

Px[_x{ = J r^rJ^{E,^,@) d cos e/Y, J p^:^+(^cm,e)dcosG. (i7) 

A± 

We note that the xf ^-iid two-body decays do not suppress the cos (j) terms, while the highest 
cos(2(^) mode in the case may be strongly suppressed due to a small polarization analyzing 
power of the decay. In particular, the coefficient becomes vanishing if two states, and ui, 
are degenerate. 

For our numerical demonstration in the xt cases, we take the masses of xt/W^ and 

yi/vii in our BCMM scenario to be 

BCMM : m± = m-± = m^± = 300 GeV and tuq = m^^ = m^^^ = 200 GeV. (18) 

Even though the i^ei/iyg are not the LSP/LKP, they decay to neutrinos and the LSP/LKP, neither 
of which is visible in the detector. 



We perform fits of the azimuthal angle distributions obtained with an integrated luminosity 

of 500 fb"^ to (1 + Aicosc/) + cos2(^)/27r. The results for the BCMM mass spectrum are 

displayed in Fig. \5\ Error bars corresponding to the 1-a uncertainty range are obtained with 

Br(xi^ — > -^^j^e) = Br(VFj^ — > i^i'n) = 0.4. First, we note that the coefficient A2 is too small 

(less than 0.5%) in magnitude to distinguish the spin-1 state from the spin-1/2 x^ state. This 

strong suppression in the case is not only due to the small analyzing power of the two-body 

decays in the BCMM scenario but also due to the cancellation of the corresponding production 

helicity amplitudes, ~ m? ±/E'^^, that is forced by electroweak gauge invariance to save the uni- 
"^i 

tarity [13] . On the contrary, the coefficient Ai in both the Xi^ and the case is sufficiently large 
so that they can clearly be distinguished from the spin-0 case. 

To conclude. The fully-correlated azimuthal-angle correlations encoding quantum interference 
between different helicity states can provide a powerful method of spin measurements at the ILC. 
We have found that, if all the particle masses are known, all the cosines of the azimuthal angle 
between two decay planes in the process ([TJ are unambiguously determined despite the inherent 
2-fold discrete ambiguity in determining the four-momenta of the decaying particles. 

Quantitatively, we have shown in a specific scenario that the spin-0 smuons can be distinguished 
from spin-1/2 KK muons or higher-spin states. However, it turned out to be difficult to establish 

''A detailed derivation of the explicit forms of the production density matrices will be presented in a separate 
publication. 
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Figure 5: Coefficients Ai and A2 as a function of Ecm for the SUSY and the UED 
production in the BCMM spectrum with an integrated luminosity of 500 fhT^ . Error bars obtained 
with Br{xi — > i^i'i) = Br{W^ ^^i^ii) = 0.4 correspond to the 1-a uncertainty range. 

the spin-1 nature of the KK Ty-boson due to the strong suppression of the highest cos 2<j) mode, 
requiring other methods such as the decay polar-angle distributions [8] and the singly-correlated 
azimuthal-angle distributions [H]. However, the latter suffers from the two-fold ambiguities in the 
reconstruction. 

Before closing, we emphasize once more that, although applied only to the SUSY and UED 
processes explicitly in this report, the proposed spin-determination method through quantum 
azimuthal-angle correlations can be used for any process with a generic event topology similar 
to that in the process ([1]) in the SM and beyond. 
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